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Abstract 

The propagation of classical gravitational waves in Bianchi Type-I universes is stud- 
ied. We find that gravitational waves in Bianchi Type-I universes are not equivalent 
to two minimally coupled massless scalar fields as it is for the Robertson- Walker 
universe. Due to its tensorial nature, the gravitational wave is much more sensitive 
to the anisotropy of the spacetime than the scalar field is and it gains an effective 
mass term. Moreover, we find a coupling between the two polarization states of the 
gravitational wave which is also not present in the Robertson- Walker universe. 
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§1. Introduction and Summary 

A Bianchi Type-I (B-I) universe, being the staight forward generalization of the flat 
Robert son- Walker (RW) universe, is one of the simplest models of an anisotropic universe. 
Unlike a RW universe which has the same scale factor for each of the three spatial direc- 
tions, a B-I universe has a different scale factor in each direction, thereby introducing an 
anisotropy to the system. It moreover has the agreeable property that near the singularity it 
behaves like a Kasner universe even in the presence of matter and consequently falls within 
the general (classical) analysis of the singularity given in jl]. And in an universe filled with 
matter for which p = 7p, 7 < 1, it has been shown || that any initial anisotropy in a 
B-I universe quickly dies away and a B-I universe will eventually evolve into a RW universe. 
Since the present day universe is surprisingly isotropic, this feature of the B-I universe makes 
it a prime candidate for studying the possible effects of an anisotropy in the early universe 
on present day observations. 

Curiously, in light of the importance of B-I cosmologies the general behavior of gravi- 
tational waves (GW) in a B-I universe has not been fully analysed. The propagation of a 
specific GW in a B-I universe has been studied before in ||. This analysis, however, was 
only done for a single wave propagating along the symmetry axis of the 2/3, 2/3, —1/3 axial 
symmetric Kasner spacetime. More recently GW in B-I universes were studied by Miedema 
and van Leeuwen within the context of a general perturbative analysis of the B-I uni- 
verse. They found, however, certain subtleties with gauge freedom and gauge fixing which 
we have not encountered. In particular, they found that an initially transverse GW will, as 
it evolves with time, become longitudinal. This we did not find and we believe that their 
result is due to incorrect gauge fixing. We shall comment further on this after our analysis. 

The purpose of this paper, therefore, is to begin the study of GWs in general B-I space- 
times. The approach we shall take follows that given in || which analyses the propagation 
of GW in the RW universe. There are, however, certain subtleties in the propagation of 
GW in B-I universes which do not occur in RW universes that makes this analysis much 
more difficult. These involve, in part, the choice of the appropriate gauge for the GW 
as well as the definition of the polarization tensors for the wave. We find that the usual 
transverse-traceless and synchronous gauge conditions that are valid in the RW universe 
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are inconsistant in the B-I universe. As it is the gauge conditions which determine the 
properties of the polarization tensor, modifications to the gauge conditions were needed to 
produce a polarization tensor which has the usual properties expected of a GW. (Indeed, if 
we take the usual transverse and synchronous gauges we would end-up with a GW which is 
not traceless.) Because, however, the usual gauge choice is no longer valid in a B-I universe, 
additional gauge dependent terms are introduced into the lagrangian which will, in general, 
cause a coupling between the two polarization states of the GW. 

The second subtlety we have encountered in defining a GW is in the definition of the 
polarization tensors themselves even after a gauge choice has been made. As usual, the 
polarization tensor can be expressed in terms of the appropriate tensor product of two 
polarization vectors which are transverse to the propagation direction of the GW. There is, 
however, a certain amount of freedom in the choice of these polarization vectors even after 
requiring that they be orthonormal to one another. Namely, one can always do a rotation 
of the polarization vectors, which can be time dependent, in the plane perpendicular to the 
propagation direction of the GW and the resultant vectors would still be a valid choice of 
polarization vectors for the GW. While this rotational freedom does not play a role in the 
isotropic RW universe, the B-I universe is anisotropic and we have found that an arbitrary 
choice of polarization vectors will cause a coupling between the two polarization vectors. 
This coupling is fictitious, however, and vanishes once an appropriate choice of polarization 
vectors has been made. In fact something quite similar happens even with GW propagating 
in RW universes if we choose polarization vectors which are not orthogonal to one another. 
This introduces a coupling between the two polarization states which vanishes under an 
appropriate rotation of polarization vectors. 

In general, we find the propagation of GW in B-I universes to be very much different 
than in RW universes. In the RW universe the two polarization states of the GW decouples 
from one another and the lagrangian for GW is equivalent to two minimally coupled massless 
scalar fields. Neither is true for GW in B-I universes. Due its tensorial nature, the GW is 
much more sensitive to the anisotropy of the B-I universe than a scalar field is and it gains 
an effective, time dependent negative masslike term. Moreover, the two polarizations of the 
GW now coupled to one another. As this coupling comes in part from the gauge dependent 
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piece of the lagrangian, its physical relevance is questionable. In particular, it is possible 
that with a more clever choice of gauge this term will disappear. After our analysis we shall 
present arguments for the physical relevance of the coupling terms. 

The rest of the paper is divided into five parts. In §2 we shall review the basic properties 
of a B-I cosmology paying particular attention to the Kasner and Zel'dovich universes. In 
§3 we shall address the question of gauge fixing and the choice of polarization vectors for 
the GW. In §4 we shall derive the equation of motion for the GW by expanding the GW in 
plane waves and writing the action for the GW in momentum space. Then in §5 we shall 
solve the evolution equations for the special case of a GW propagating along an asymmetry 
axis in a Kasner universe or a Zel'dovich universe. This is the only case in which closed form 
solutions of the equations of motion can be found. Fortunately, they are the most physically 
relevant. Concluding remarks can then be found in §6. Finally, in the appendix we shall 
repeat the analysis in §4 using a different gauge choice for the gravitational wave and show 
that coupling terms similar to those found in §4 appear even in this gauge choice. 
§2. Review of B-I Cosmology 

In this section we shall present a brief review of B-I cosmologies in the absence of per- 
turbations. In doing so we shall follow the notational and sign conventions found in ||. In 
particular, one may always choose coordinates for the B-I spacetime such that the metric 
has the form 
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ds 2 = dt 2 -Y^a-^idx 1 ) 2 , (1) 

i=l 

where as usual Greek indices will run over the four spacetime directions while Roman indices 
will run only over the spatial directions. Although we shall use the summation convention 
for Greek indices, for clarity we shall not use it for Roman indices. 

It is convenient to define a pseudo-conformal time variable r\ through 

^ = -L (2) 
dt a( V ) ' 1 ' 

where a = (010203) is the geometric average of the three scale factors along each direction. 
Then the only non-vanishing components of the Riemann tensor are 



while for i j, 



where C = a 2 (rj), Cj = a 2 Arj), 



1 didjCiCj , . 



C C 1 3 

and the prime denotes derivative with respect to rj. Einstein's equations in the presence of 
matter are then 

-£>' + eg = -4ttC(p + 3p) , 



where 



d'j + djD = 8nC(p - p) , (6) 



Q^^Eidi-dj) 2 , (7) 

IZ i>3 



and 



( 1 da\ 2 87ra 2 ^ ,„ s 
(a^j =— " + <? - < 8) 

This, aside from the term dependent on Q, is identical to the RW case. Q is thus the physical 

measure of the anisotropy of the spacetime and is called the anisotropy factor. Taking then 

the difference between two different directions in eq. (|6]), we find that 

di-dj = (d° -d])/a 2 , (9) 

where d® is the value of di at some initial time r] for all p and p and we have chosen the 
overall scale such that a(r] ) = 1. Consequently, Q = Qo/a 4 where Q = Q(r] ). 

Eq. (^|) has been solved in the presence of various types of matter || . We shall, however, 
only be concerned with two special cases in this paper: the Kasner universe, which is free 
of matter, and the Zel'dovich universe where p = p. In both cases, 

a 2 = y/y , a 2 = a 2 (y / 'y ) 3p * , (10) 



where y = 1r\^[Q~§ for the Kasner universe while y = 2r]y8irp /3 + Q$ for the Zel'dovich 
universe, po is the energy density of the universe at rjo while a^j = aj(0) are the initial scale 
factors in the various directions which, without loss of generality, can be set to unity. 



Pj are parameters of the B-I spacetime which measure the relative anisotropy between 
any two asymmetry axis. As they must satisfy the constraints 

1=X>, ! = !>?> (11) 

i=i -i=i 

out of the three parameters, only one is arbitrary. Since eq. (|11]) describes the intersection 
of a sphere with a plane in the parameter space (pi,P2,P3), we may parametrize the allowed 
values of pj by an angle on the unit circle 0. One particular choice of parametrization is 

Pi = — (l + cos 9 + v^sin^ , 
P2 = - (l + cos 6 — \/3sin#) ( 

p 3 = ^(l-2cos6) , (12) 

Although a priori 9 ranges over the unit circle, note that the labeling of each pj is quite 
arbitrary. Thus the unit circle can be divided into six equal parts each of which span 60°, 
and the choice of pj is unique within each section separately. Notice that when 6 = 0, 
Pi — Vi = 2/3 while p% = —1/3. This is the spacetime considered by Hu in ||. When 
9 = 7r/3, on the other hand, p\ = 1 while p<i = p% = 0. It can be shown that this spacetime 
is equivalent to the Minkowski spacetime up to a coordinate transformation. 
§3. Gauge Conditions and Polarization States 

In this section we address the problem of gauge fixing for the GW and the subsequent 
definition of polarization states. As usual, we consider the GW as a perturbation off the 
background metric by writing g^ u = g® u + h^, where g Q is the unperturbed B-I metric 
given in eq. (0) while is the perturbation which has a gauge freedom to be taken care 
of. The usual gauge choice, which works for both flat Minkowski and RW spacetimes, is 

V"V = ° > « M V = ° > &£ = 0, (13) 

and are called the transverse, synchronous and traceless conditions, respectively, is a 
timelike Killing vector which, without lose of generality, can be taken to be (1,0,0,0). In 
the RW spacetime, the transverse and synchronous condition actually implies that the GW 
is traceless and this condition is redundant. Unfortunately, for a B-I universe this is no 
longer true. 
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Let us express hfy as a plane wave 

h^^Ml)^. (14) 

where q ■ x = J2j Qj^ an d s labels the polarization state. Because g Q is a function of 77 
only, we can always make this choice. Our convention is that the qj are independent of 77 
while qi = g^qj. zuff is the polarization tensor for the GW and must, by definition, be 
independent of the amplitude h^ s {rj) of the GW. The usual transversality condition then 
gives 

= <r + ™ {s) od» log h* + * E ffM? + E - E , (is) 

as well as 

= gTd^ + wMfa log ^ + i E + . (16) 

j a 

From this, we see that for the polarization tensor to be independent of h$ s , we must choose 
the synchronous gauge. After doing so, the transversality condition reduces to 

o = E-^ > o = e^M?- ( 17 ) 

3 a i i 
and becomes the standard transverse-traceless condition in a RW universe. For a B-I uni- 
verse, on the other hand, we see that the GW is no longer traceless if we use the usual 
transverse and synchonous gauges. The usual gauge choices eq. (|l^) are no longer self- 
consist ant in a B-I universe and must be modified. 

There are two straightforward modifications that we can make. The first is to use the 
usual transverse and synchronous gauges and live with a GW which is not traceless. This, 
however, introduces a great deal of complexity to the problem and we shall not do so 
(see appendix). The second is to require that the GW be traceless, but to modify the 
transversality condition. Namely, to choose 

= V^ , o = v > = /#, (18) 
which gives the following constraints on the polarization tensor 

o = E<?M? , o = tu« , o = 5><->;. (19) 

3 3 
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Notice that these equations determining gjW have precisely the same form as the those in 
either a RW or Minkowski spacetime which is the reason we shall work primarily in this 
gauge. 

(s) 

Instead of the polarization tensor w\j , we shall find it conceptually easier to work with 
the polarization vectors with 

„(+) _ J2) (2) (x) _ (1) (2) (2) (1) r?m 

as long as 

o = E^'4 s) , E4 S) ^' )J = -^'- (2i) 

Implicit in this definition is the added requirement that the two polarization states be "or- 
thonormal" to one another: wfym^ 1 *" = 25 s s> . Notice, however, that this orthonormality 
condition is not required by the gauge choice eq. ([19]) but is rather made separately as a 
minimal requirement for the two polarization states to decouple from one another. While 
such a choice is sufficient for a RW or Minkowski spacetime, we shall see in the next section 
that additional requirements are needed for a B-I spacetime. 

Eq. (|2lD does not determine ef^ uniquely, however. There is still a rotational freedom in 
the plane perpendicular to qj. Namely, we can rotate 

"(!) (!) i I (2) • , 

€j = £j COS (p + €j Sin , 

ef = -ef sin + ef cos , (22) 

and 6j will still satisfy eq. (|2ip . In particular, may also be ^-dependent. Moreover, 
once a specific choice of polarization vectors has been made rotational symmetry will be 
broken. While this breaking of rotational symmetry is of no consequence in the RW universe, 
which is isotropic to begin with, it is of great consequence in the anisotropic B-I universe. 
In particular, any specific choice of polarization vectors will tend to break the expected 
exchange symmetry 1-2-3 among the labels of the axis of the spacetime. We shall also find 
that the equations of motion for the GW will be greatly simplified if we make a special 
choice of polarization vectors. For these reasons we shall work as much as possible with a 
general set of polarization vectors and shall delay making a specific choice of polarization 
vectors until forced to. 
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To find a specific representation of we first define a local coordinate system using the 
vierbeins e* , 

e^w = -S ij . (23) 

Then, taking = (0,ai,0, 0), = (0,0, a2,0), = (0,0, 0,03), one specific choice of 
polarization vectors is 

(1) 9l<?3 1 . 92<?3 2 (To I -2 3 



i 2 ) - 92 ^1 gl ^2 



, 9l + 92 V9i + 92 

where 

q 3 = - jMtl . (25) 
Notice in particular that ef\—q) = (—l) s+1 €j S \q) so that w^{— q) = m^){q) while 

The polarization tensors defined in eq. (|24|) are for linearly polarized GW. We shall also 
consider circularly polarized GW defined through 

_ 1 r„(+) , ,-„(xn _(*) _ 1 /_(+) ■ (x)\ , 9fn 

with w^ L >[— q) = zu^(q) = w^ L \q) and vj^ R \—q) = vj( L '(q) = w^ R \q) and the bar denotes 
complex conjugation. We can, of course, also define circularly polarized polarization vectors 
in precisely the same way 

£ f> = -L( £ « + i£ f>) , £ <*) = _L (£ a)_ i£f)! (27) 

and there is a particularly simple relationship between w and e, 



EC 



(L) " , = el R \ (28) 



Finally, we note that under a rotation of polarization vectors as in eq. (|22"D , e*- = e l< ^- ', 



54. The Action 
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In this section we shall derive the equation of motion for the GW. The approach we shall 
follow was first used by Ford and Parker in their analysis of GW in RW universes and 
involves the expansion of the GW in plane waves. To lowest order the action for the GW is 




which differs from the one used in |J for GW in RW spacetimes by the addition of a 
kinetic term. This is because our gauge choice eq. (|Tj|) is different from the one they used 
as in eq. flllf ). We should also mention that even if we did use the usual transverse and 
synchronous gauges, our action will have additional terms dependent on h)t which no longer 
vanishes in this gauge for a B-I universe. Indeed, any choice of gauge for the GW in a B-I 
universe will introduce additional terms to the lagrangian. The form that these additional 
terms take is dependent on the choice of gauge one makes, however, and at first glance the 
physical relevance of these terms, since they are seemingly gauge dependent, is questionable. 
We shall see, however, that although the explicit form of these terms changes with the choice 
of gauge, they all have the same root physical cause. Namely, they are due to the fact that 
in a B-I universe one cannot define a polarization state for the GW which will not change 
with time. Consequently, these terms are physically relevant. 

It is straightforward to show that due to Einstein's equations, the second and third terms 
in eq. (^9) cancel one another. We are then left with three terms to consider 

Ik = \J aWxVX^h* , 
I R = -^fa 4 d 4 xR^ au h^, 

h = -\\ a 4 d 4 xV p h^V a h« , (30) 
where I = 1% + Ir + I g . We next expand h pu in plane waves, 

V = £ E (3i) 

q s=+,X 

where h s (q, x) = h^ s e Lq ' x with the reality condition given by 

h+(q,x) =h+(-q,x) , h x (q, x) = -h x (-q, x) . (32) 
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Turning our attention first to Ir, we use eq. ([IT]) and find that 

J R = ~\^1Y^ J a 4 d 4 xT ssl (e)h s ,(q,x)h s (q,x) , (33) 



q s 

where 

T ss ,(e) = Rf^w^wff (34) 
In terms of polarization vectors, 

T ++ (e) = T xx (e) = -2 £ i^ e W e Mef e< 2 > , (35) 

ijkl 

while 

T +x (e)=T x + (e) = 0, (36) 
due to the anti-symmetry properties of the curvature tensor. Next, using eq (f|), 

where it is important to note that we have not as yet chosen a specific polarization vector 

e (s) 

We next consider the gauge term which becomes 

/ S = ~EE/ aW^'W^ V^V<A + 2s7£V a t*«' ltt V'7U, + 

(f ss' \ 

V?w$V a w^ m h s ,hX , (38) 



for plane waves. The first two terms involve the term zu^ ^ a V a h s (q,r]) which vanishes 
identically due to eq. flT9|). The last term does not, however, and if we define 



= -V^V-tuW; (39) 

we find that 

* j 
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(i)' 



1 /^a- 

E— ro 



> i ^ -- (x) 2 



j J 



^eNe^ii (40) 



o I / j I \ \ Z j 7 

GT V ; Cli J \ j dj 



while N+ x = N x + and we have used V^ej^ = 0. Notice, in particular, that in general 
N+ x 7^ and this gauge term introduces a coupling between the two polarization states. 
In addition, N+ + ^ N x x . Notice also that each of the N 2 are related in a simple manner 
to the trace of the polarization vectors which would vanish under the usual choice of gauge 
eq. (0). 

Finally, we turn our attention to the kinetic piece, the most difficult one to work with. 
Proceeding exactly as before, we find that 



Ik = \T. J a i d A x\^V p h s VPh s -Y,Mirh sl h s -2Y,D^h s ,V^h\ , 

q \ s ss' ss' ) 



(41) 



where 



D M = -iwWJV^f , Ml, = -\V P ^}V^ V . (42) 
This can be reduced to manageable form by using polarization vectors: 

M 2 ++ = M 2 X = £ VpeWV'eM" _ 2 (e^V p e^)(e^e^) , (43) 

s 

where we have used ej^eW = — 5 SS ', while M +x = 0. Similarly, 

D = D = e (1) V e (1) " + e (2) V e (2) ° (44) 
which vanishes identically while 

D — —D = e«V e (2) ° - e (2) V e (1) ° (45) 
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This term will not vanish in general, however, and we see that the kinetic piece of the 
lagrangian also introduces a coupling between the two polarizations. Notice also that the 
only non- vanishing term in D^ +x is D 0+x = D +x . 

Combining these three pieces together, we finally arrive at 

-- 2 D x+ h + d ^- - - 2 D +X h x ^± + % (Mx + h x h+) } , (46) 
a 1 ar] cr drj 2 v ' J 

where 

m^ e ) = a 2 (M 2 + + T ++ - iV s 2 s ) . (47) 

Once again it is important to note that we have not, as yet, made any specific choice of 
polarization vectors. 

At this point we can see the differences between the propagation of GW in a B-I uni- 
verse versus a RW universe. Ford and Parker has shown that in a RW universe the two 
polarizations of the GW decouple from one another and each, separately, is equivalent to a 
minimally coupled massless scalar field. Neither is true for the GW propagating in a B-I 
universe, however. Here the GW picks up an effective mass term due to it being a spin-2 
field and is much more sensitive to the anisotropy than a scalar field is. What is even more 
surprising is the apparent coupling which is present between the two polarization states of 
the GW. 

The physical relevance of this coupling term is somewhat questionable at this point, 
however. They could have arisen from an inappropriate choice of polarization vectors or, 
since the second such term comes from the gauge piece of the action, because we made an 
inconvenient gauge choice for the GW. For example even in a RW universe a coupling term 
between the polarizations would appear if we had chosen polarization tensors which were 
not orthogonal to one another; an "inappropriate" choice of polarization tensors. To begin 
to separate the truly physical effects of the anisotropy on the GW from the arbitrariness in 
defining the GW, we have to develope a better understanding of the effects of the anisotropy 
on the GW. This is best done by considering the behavior of the polarization vectors instead 
of the full tensors. 
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Consider the triad (jjj,6j ,6j) which form a local orthonormal coordinate system. Be- 
cause we are in the B-I universe, these three vectors are not fixed, but rather changes with rj. 
Since they are constrained to lie on the unit sphere, their motion consists of two rotations; 
one in the plane perpendicular to q which is spanned by ep and ef \ the other parallel to 
q. Let us consider the two rotations separately. Rotations parallel to q occur because the 
direction of propagation of the GW, q~j, is always changing with time (see eq. (250). This is 



because the medium through which the GW is propagating, the B-I universe, is anisotropic 
and always changing with rj. In particular, notice that in a RW universe, where such cou- 
pling between polarization states is not present, cjj does not change with time and there 
are no rotations along this direction. Rotations in this direction are therefore caused by a 
physical effect of the anisotropy on the GW. 

Rotations in the e^-e*- 2 -' plane, however, are not physical. They can occur even when 

(s) 

q~j does not change directions as in a RW universe. Remember that the definition of el- 
is somewhat arbitrary. One still has the freedom to do a rotation as in eq. (^) of the 
polarization vectors in this plane even if this rotation is //-dependent. Returning to eq. fl45|) 



we see that D +x = 2e^Vo^ 2 ^ — ~ 2e^ Voe 1 ' 1 ^ depends explicitly on the velocity Vo£j of 
the polarization vectors in this plane. Consequently, the presence of the D +x coupling term 
is due solely to the rotation of the polarization vectors in the plane perpendicular to qj. A 
-D+x 7^ therefore only means that we have not chosen the "correct" polarization vectors: 
one in which the polarization vectors do not rotate about qj. 

To demonstrate that such a choice of polarization vectors always exists, we preform the 
//-dependent rotation of the polarization vectors given in eq. (0). Under this rotation, we 
find that 

£> +x (e) = 20' + D +x (e), (48) 

while 

/4(e) = /^(e) - 4(0') 2 - 80'^TefVoe^ - 

3 

a 2 [N 2 ++ (e) (cos 2 20 - 1) + N 2 + x (e) sin 40 + iV 2 x (e) sin 2 2 
ml (e) = ml (e) - 4(0') 2 - 80' £ ef V e^' - 
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a 2 {^ + (e)sm 2 20-^ x ( e )sin40 + ^ x ( e )(cos 2 20- 1)} , 
N 2 + x (e) = \ (Nl x (e) - N 2 ++ (e)) sin 40 + N 2 + x ( e ) cos 40 , 



(49) 



Clearly, one needs only to choose for any given e - 



' +x (e)dr]' + 0o , 



(50) 



and the coupling term proportional to D +X (e) vanishes. 0o is an arbitrary, constant angle 
and its presence means that we still have a degree of freedom to choose our polarization 



Notice also that eq. ([51]) does not have a symmetry under exchange of 1-2-3 which is a reflec- 
tion of the fact that once a specific choice of polarization vectors has been made rotational 

(s) 

symmetry is broken. Other choices of e) will only result in different forms for 0, but D +x 
will still vanish. 

While D +x vanishes after an appropriate choice of polarization vectors is made, this still 
leaves the coupling term proportional to N+ x . This term came from the gauge dependent 
term of the original action, however, and is different for different choices of gauge. It is 
possible that under other, more clever gauge choices this term will no longer be present. Its 
physical relevance is therefore questionable at this point. There are two ways to address 
this problem. The first is to redo this calculation with a different gauge choice and see if 
the coupling terms are still present. This we have done in the appendix. This method, 
however, has the shortcoming in that one can never be sure that if there is still some other 
gauge choice for which the coupling terms vanish identically. The second way is to develope 
a physical argument for or against the relevance of the physical terms. The final conclusion 
of this argument should therefore be valid with any gauge choice. This is the approach we 
shall now take. 

Consider the function 



vectors. For the specific choice of polarization vectors given in eq. (p4"D, 




(51) 




(52) 



15 



where 771 is some fixed time and e^ J (971) = g 30 {r]\) e^j(r]i). Clearly, f ss '(r]i) = —5 SS / and thus 
f ss r is a measure of how the orthonormality of the polarization vector changes with time. 
Next, considering f ss t as a function of rj only, we expand €j(r]) in a Taylor series in rj about 
rji and find that to first order, 



fss'(v) = Sss> + \J2 e 



(s)3 de j 



tn - Vi) + ■ ■ ■ ■ (53) 



dr. 

Since rji is arbitrary, we see that the infinitesimal change in the normality condition is 
measured by 

^ c (l)3 ae 3 = ^ a j c (l)3 c (l) (2)3 a6 j = a 3 (2)3 (2) ^ 

y dv 3 ' y ^ i a i J 



The difference of these two terms is just ^—a?N+ + = J2^ + ^j(^j/ a j while their sum is just 
the T j a 'j/ a j defined in eq. ( |A2| ) of the appendix. For the choice of gauge eq. flTg| ) it plays 
no role. The infinitesimal change to the orthogonality condition is measured by the linear 
combinations 

£ i^'^h + ^'^h I = 2£- e(1) 'ef , (55) 



which is just J—a 2 N' : 



r 2 

X X 




(56) 

We thus see that each of the terms which may cause a coupling between the two polariza- 
tion states, N+ + , N xx , N^_ x , and D +x , have their roots in the time rate of change of the 
polarization vectors. Since the polarization vectors change with time due to the anisotropic 
rate of expansion of a B-I universe, we therefore conclude that the presence of these coupling 
terms is a natural consequence of the physical properties of the spacetime. In particular, 
we see that a coupling between the two polarization states will be present so long as we are 
not able to consistently define a polarization state which is valid at all times. They will be 
present in some form no matter what gauge conditions one chooses and the coupling terms 
are physically relevant. 
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To further illustrate the connection between the rate of change in the direction of the 
polarization vectors and the coupling term, let us look at the conditions under which theses 



coupling terms will vanish. For the special choice of polarization vectors given in eq. (24) 
q( + qi \ai a 2 J 

N 2 = 2 Q1Q2Q3 (<A _ 4\ 

+x a? q\ + q% \a\ a 2 J 

2 / a 2 a 3 1,-2/ fl 3 



\a 2 a 3 / \a 3 «i / <?i + 92 









a 2 J _ 



(57) 



We caution the reader that while D +x is invariant under constant rotations of the polariza- 
tion vectors, N+ x is not and our arguments are valid only within this choice e- . We see that 
all coupling between the two linear polarizations will vanish if either ai = a 2 or if qj = 
along some direction j. Although the first condition seems to break the 1-2-3 exchange 
symmetry in the labeling of the anisotropy axis, this simply means that we were fortunate 
enough to choose the "canonical" polarization vectors for the GW which decouples the two 
polarizations. Any other choice of polarization vectors will produce an N+ x ^ 0, but can 
always be made to vanish with a constant rotation of the polarization vectors. D +x , on the 
other hand, vanishes for all such choices of polarization vectors. 

Referring to eq. (|2~4| ) we see that when a\ = a 2 , tj lies in the 1-2 plane and, more 
importantly, its direction does not change with time; Vo^ 2 '' is once again parallel to e^ 2 -*. 
Thus any direction we choose for this polarization will always lie in this direction and will 
not change with time. 

Suppose, now that a\ ^ a 2 , but g 3 = 0. Then all the coupling terms between the two 
polarization states vanish. In this case qj lies in the 1-2 plane but more importantly we may 
always choose one of the polarization vectors to lie in the 3-direction perpendicular to this 
plane. Once again we see that the direction of this polarization vector does not change and 
we can once again have a consistant definition of polarization vectors for all time. 

We now see explicitly that the coupling terms vanish as long as the direction of one of 
the polarization vectors does not change with time. This also can be seen in the definition of 
N+ x which depends on the directional derivative of the polarization vectors. Physically, it 
means that when this happens we can consistantly define the direction of at least one of the 
polarization vectors at all times. The second polarization can be found by taking the cross 
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product of this polarization vector with the direction of propagation of the GW. Although 
this can be done in certain special cases, it is not true in general and there will in general 
be a coupling between polarization states of the GW no matter which gauge one picks. 
Finally, we note that for circularly polarized GW 

/ = £ / aWx I £ 1 (v p h s V p h s - ^fih s h s ) + ^h R h L + ^h R h L \ (58) 

q [s=R,L \ CI J CI CI J 

where h R = (h + — ih x ) / \^2, h L = [h + +ih x )/ y/2, 

(3(e) = e^(3{e) , (59) 



with 



P(e) = \ [e«V^ 2) + efV^ + i(e«V^ - 6«V%«)1 , (60) 



and (p defined through eq. (p0|). We have thus eliminated the D +x with this choice of 
polarization vectors. Notice that for circularly polarized GW the mass terms for both 
polarizations are equal 



m 2 (e) = ^[m+(e) + m 2 (e)] 



= - £ (■* - ?Y # - E - + 4v>Lm + *!>.(«» , (en 

cyclic \ a 3 a k) j>k\ a 3 a * / 2 

but we pay for this symmetry through the addition of a complex coupling term between 
the two polarization states. This coupling term will only vanish when both the real and 
imaginary parts of /3(e) vanish. If, however, they vanish for one choice of polarization 
vectors, they will vanish for any other choice from eq. (|59"|). Consequently, we may use the 
polarization vectors given in eq. (0), and find that 

= 6« V-4 2) + *fV»$> = -j±k= fe-^), (62) 

as well as 

= cWVe^ - tf'Veg = -g 2 f^ 2 - + g - 4) + ^ 
^ \ci2 a 3 J \a 3 a x ) qf + q£ \ai a 2 J 

(63) 
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must both vanish separately. This can only happen when the GW propagates along an 
asymmetry axis and the spacetime is axially symmetric in the plane perpendicular to this 
axis. For example, the coupling term will vanish if q 3 = — 1 and = a 2 - This is a much 
more stringent condition than for linearly polarized GW and is because although the two 
linearly polarized GW may decouple from one another, they will still have different mass 
terms and this introduces an additional coupling. Thus we have the peculiar situation that 
while the linear polarized states may decoupled from one another, the circularly polarized 
states need not. This once again underscores the fact that in a general B-I universe one 
cannot consistently define a polarization state at all times. 
§5. Solutions of Equation of Motion 

The evolution equation for GW in an arbitrary B-I spacetime cannot be solved in terms of 
known functions in general, especially when there is a coupling between the two polarization 
states. Consequently, we shall only consider GW propagating in either the Kasner or the 
Zel'dovich spacetimes where solutions can be found in certain special cases. 

From the action eq. (46) the evolution equation for a linearly polarized GW propagating 



in a Kasner or Zel'dovich spacetime is 



dy 2 a dy dy 
d 2 h r 7* _ 1 da dhr? x I " 

,< 1 




(64) 

where we have not chosen a D +x = and we have taken y — 1 for convenience. In the above, 
the tilde denotes the fact that we have scaled the corresponding quantity by either 4Q or 
4(87rp /3 + Qq) depending on whether the spacetime is a Kasner or a Zel'dovich universe 
to make the resulting quantity dimensionless. For example, q 2 = q 2 /4/Q for the Kasner 
universe while q 2 = g|/4/(87rp /3 + Q ) for the Zel'dovich universe. We see that due to the 
coupling term between the two polarizations, even if the GW is initially plus-polarized, a 
cross polarization will be generated. Unfortunately, even in the case of the Kasner universe 
eq. ([64]) cannot be solved generally. We shall instead have to look at limiting conditions. 

Suppose, for convenience, that pi > p 2 > p%. Let us now consider a general q for which 
each of the components of q do not vanish. If we then take the limit y —>■ we find that 
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qi = —1 while q 2 = qs = 0. With the choice of polarization vectors given in eq. 



a' ~'^ 2 



ml = , m 2 x = -p-^ , (65) 

\a 2 a 3 j 

where the prime now denotes the derivative with respect to y. N +x = D +x = and there 
is no coupling between the two polarizations. Similarly, when y — > 00, qi = q 2 = while 
q 3 = — 1. Now 

ml = , m 2 x =-f^-^V. (66) 

In these two limits the behavior of the GW simplifies dramatically and can be solved in 
closed form. More importantly, these two limits, which corresponds to the behavior of the 
GW near the initial singularity as well as the large time behavior of the GW, are the most 
physically significant. 

It is therefore sufficient to consider a GW propagating along one of the symmetry axis; 
say in the /th direction. Then in eq. (|6^) and h^ x decouple into two equations 

_ d 2 h ?+ lda dhfr 1-3p 2 
dy z a dy dy 

= fhf + + L-*^ _ (<k _ £f*VX , Vx , ( 67) 

dy z ady dy \aj 

where I ^ j ^ k. Notice that in general the behavior of will be different from that 
of hfx due to the additional masslike term in its equation of motion. It is only when we 
are propagating along the asymmetry axis of an axially symmetric B-I universe will the two 
polarizations have the same behavior. Since the solution for h$ + can be obtain from that of 
hqx by setting aj = a& in the mass term, we shall concentrate our attention on the solution 
of h$ x . 

It is straightforward to see that the solution to eq. (|67|) are Bessel functions J u and Y v . 



Using the parametrization of the pj given in eq. (12), we find that 



h^(q,x)=e^ xl \A^J u+ 



J->x ± y+ 




1 - \ cos 9 - ^ sin e 
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{/ 7i „,l+cos 9 \ 
^ (m) + 

,l+cos0\ 1 



1 + cos6>^ 
where 

13 cos 9 =F a/3 sin 9 



Bfy, q ^— ] ) . (68) 



2-cos^T V3sin9 ' 
/- I sin ^ I 

"3 = ^3-i- — , 69 

1 + cos 9 

while and are integration constants. The superscript signifies that this is a 
GW which propagates along the jth direction. The solutions are obtained from the 
above by replacing the Bessel functions of various orders with Bessel functions of order zero 
but with the same arguments. 

Of particular interest is the small and large y limits to eq. (|6S|). The large y limit deter- 
mines whether or not the GW (which is a perturbation of the metric) increases or decreases 
in magnitude. It thereby establishes whether or not the Kasner and Zel'dovich universe 
are stable under small perturbations. The small y limit, on the other hand, determines the 
behavior of the GW near the initial singularity. To do so, however, we shall have to find a 
way of comparing matrices. We therefore define the norm 

\\T^ U \\ = max \T^ V \ , (70) 

where the max is taken over all components of the tensor. This norm has the advantage of 
having all the nice properties associated with a norm ||. 
Taking first the y — > oo limit, we find that 



\ 



2 - cos 9 = V3smg (2 _ cos g_y3 Sin g )/4 



4 1} cos^ 



(2-cos6»-v / 3sin6»)/2 

2 — cos 9 — a/3 sin 9 
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W + + 2j +5 >< Sm 2-cos^- v / 3sin^-2r + 2 



(71) 



while 



q3 yl+cos8 n 



hV&x) ~ / 2(i + co S g) (1+cose)/2 r (3) cog _ ^ + 

R (3) . ( 2q,y^ e 7T / , 1\\ \ 

Bx sm (TT^-2r + 2jJ|' (72) 

with /i^ being obtained from by taking — > —9 and replacing (1) with (2) every- 
where. It is then straightforward to show that the amplitude of ||ftfe;||/||<7|«/|| decreases with 
increasing y. Consequently, we see that the perturbation of the metric decreases with 
increasing y. 

Taking now the small y limit, we find that 

h (s) „ y -V-^e ) (?3) 

where once again is obtained from by taking 9 — ► —9. Eq. (|73|) holds as long as 
the degree of the divergence in does not vanish. When it does, that will have a 
logarithmic divergence. We therefore find that as long as P2 ^ P3, 11^^11/11^11 ~ 1 from 
GW propagating along the 1 or 3 direction while ~ for a GW propagating 

along the 2 direction as y — > 0. Similarly, H/i^ll/llflwII ~ l°gl/ f° r a GW propagating along 
the 1 or 2 directions and ||fr^||/||(7/xi/|| — > for a GW propagating along the 3 direction. We 
thus see that along the 1 and 2 directions the perturbation of the metric becomes unboundly 
large near the singularity. 

For the special case of P2 = P3, we find that ll^^ll/II^H ~ 1 for GWs propagating 
along the 2 or 3 directions while W/Wg^uW ~ logy for a GW propagating along the 1 
direction. Similarly, we find that H/i^ll/lls^H ~ logy for a GW propagating along the 1 
or 2 direction, while \\h^) W/Wg^W — ► for a GW propagating along the 3 direction. Thus, 
in this special case the perturbation of the metric becomes unboundly large along all three 
directions. The case considered by Hu, which corresponds to 9 = 0, belongs to this catagory 
and the results we have obtained agree with his. 
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The case where 9 = n/3 is quite special since this corresponds to p 2 = P3 = while 
Pi — 1. Supposedly, this case also falls within the above analysis and we would still obtains 
the same small y behavior. When p = this is the Rindler spacetime, however, which is 
known to be equivalent to Minkowski space and supposedly the propagation of GW in a 
Minkowski spacetime does not have a singularity. To address this problem, we now perform 
a more detailed analysis of the solution in this special case. 

Going back to the equation of motion, we see that when 9 = ir/3, 

h%\ qi ,x)=v£\ qi )hU 



= {4el ~ elel) [A { ? J Q (q 2 t) + B^Y (q 2 t) 



h^{q 3 ,x) = w^{q 3 )hf 



if A ~ ejej) [Af J (q 3 t) + B®Y (q 3 t) 



o l Q3 x 



(74) 



where we have put back in the polarization tensors explicitly, t = 2i] /3 and we have used 
y = (t/t ) 2 / 3 for the Rindler spacetime. Similarly, 

h%\q 1 ,x)=w$(q 1 )hy 

= (e 2 e 3 + e 2 e 3 ) e iqi< - xl+tologt/to) + e iqi ^ xl ~ tolost/to) ) 

= - (e^el + e l u el) {A { ? J, (q 2 t) + B^Y X (q 2 t) 
h^{q 3 ,x)=w^{q 3 )hf 



- (iel + elel) [A^J, (q 3 t) + B^Y, (q 3 t) 



and we can see explicitly the singularity in the solutions when t — > 0. 

Let us now do a coordinate transformation into Minkowski spacetime. When 9 
the metric is 

ds 2 = ( dt y _ (T) 2 (^1)2 _ (dx 2 ) 2 - (dx 3 ) 2 . 
To map this into the Minkowski spacetime, we make the coordinate transformation 



(75) 



= 7T/3, 



(76) 
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T = t cosh^V^o) , X 1 = t sinh(x 1 /to) , X 2 = x\ X 3 = x 3 . (77) 

Notice, however, that because h^ u is a tensor, it transforms as 

dx a dx 13 

^ = dX^dX^ P ' ^ ' 

where m^V> is the transformed GW. Note, however, that dx a /dX^ is block diagonal and 
mixes the 0-1 components of the GW while leaving the 2-3 components alone. For GW 
propagation in the 2 or 3 directions, this means that the transformed m^-oi/ 7^ 0, although 
q 3 M hij = still and it will still be traceless. Consequently, for GW propagating in these 
two directions, the transformed GW will not satisfy the usual synchronous gauge and is not 
what one usually calls a GW in a Minkowski spacetime. 

For GW propagating along the 1-axis, however, the situation is a little different. Notice 
first that the polarization matrices will not be changed by the transformation and that they 
are, in fact, precisely the polarization matrices for a GW propagating in this direction in a 
Minkowski spacetime. Next, using the coordinate transeformation eq. (|77|) , we see that 

A +r^S +5?) r^) }' (79) 

and the equation for the x mode has the same form. Notice that eq. ( |T9"D is a solution of 
the wave equation in Minkowski spacetime, as a GW in a Minkowski spacetime should. In 
this sense it is a GW in the Minkowski spacetime, but it is not plane wave and satisfies 
different boundary conditions. Consequently, we see that after a coordinate transformation 
into a Minkowski spacetime, the GW in the Rindler spacetime will not be transformed into 
what we usually refer to as a GW in the Minkowski spacetime. 
§6. Concluding Remarks 

To conclude, we have begun the analysis of the propagation of GW in B-I universes using 
the method developed by Ford and Parker. We find the behavior of a GW in a B-I universe 
to be very much different than a GW in a RW universe. The two polarization states are not 
equivalent to two minimally coupled, massless scalar fields. Rather, each polarization state 
gains what is effectively a time dependent mass term due to the tensorial nature of the GW. 
Namely, the GW is a spin-2 particle instead of a scalar particle and is consequently much 
more sensitive to any anisotropy in the universe. We have also found the two polarization 
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states to be coupled to one another. This coupling term depends explicitly on the gauge one 
picks for the GW and varies as one varies the choice of gauge. No matter what gauge choice 
one makes, however, a coupling between the two polarization states will always be present 
and the coupling is not a gauge artifact. 

The reason for this is fairly straightforward. We can consider the propagation of a GW 
in a B-I universe as the propagation of a wave through an anisotropic medium. Although the 
polarization and direction of propagation of this GW is initially arbitrary, as it propagates 
through the anisotropic medium the medium itself will tend to change the direction of 
propagation, and thus the direction of polarization, of the wave. This can be seen explicitly 
in a Kasner universe where the anisotropic medium gradually forces the wave to propagate 
only along one of the asymmetry axis. The anisotropy of the medium is constantly changing 
with time, however, and thus the polarization state that the medium wishes the GW to adopt 
also constantly changes. This is done through the coupling term between the polarization 
states and is caused by changing of the polarization vectors with time. 

We would therefore expect a coupling between the two modes of the GW to be present in 
any anisotropic expanding universe. Indeed, such a coupling has also been found by Ezawa 
and Soda who analysed the effects of the topology on the propagation of GWs. They 
considered GWs propagating on plain symmetric spacetimes with two of the spacetime 
directions compactified into a torus and also found a coupling between the plus and cross 
modes of the GW. 

Hu considered a B-I spacetime for which pi = p 2 = 2/3 while p% = —1/3 (we have changed 
slightly the notation used in |J). Moreover, he only considered a GW propagating along 
the cylindrical axis. We have seen, however, that the behavior of GWs in this spacetime is 
quite special and atypical. In fact, it is only for a GW propagating in this manner in this 
spacetime that the two linear as well as circular polarizations decouple from one another. 
Propagation along any other direction in the spacetime will introduce a coupling between 
the two polarization states. 

Miedema and van Leeuwen have also analysed the propagation of GW in B-I universes 
within a general analysis of perturbations of a B-I universe. Looking at the GW component 
of the perturbations, they found that while a GW can be defined as being transversal at any 
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given time, as the wave evolves in time non-transversal components of the wave will begin 
to appear. They thereby conclude that GW in B-I universes are in general non-transversal. 
The longitudinal components of the wave, however, were considered to have no physical 
meaning since they can be gauged away at any given time. We did not encounter any such 
subtleties in our analysis. We have instead found that with the correct gauge choice the GW 
will always be transverse. 

The difference between our two results can be found in the gauge conditions that we 
have taken. Miedema and van Leeuwen has chosen to work in the synchronous and traceless 
gauges (see eqs. (38) and (108) of ||). They then attempted to enforce the usual transverse 
condition on the GW (eq. (114) of ||]). As we have shown, these three gauge conditions 
are incompatible and one of them must be modified. By imposing the synchronous and 
traceless conditions on the evolution equations for the GW (eq. (135) of Q), they cannot 
then choose the usual transverse condition. This inconsistancy may be the root cause for 
the generation of the longitudinal components as the wave propagates that they observed. 
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APPENDIX A 



In this appendix we shall repeat the analysis found in §4 using the usual transverse- 
synchronous gauge in order to show that coupling terms between the polarization states are 
still present in this gauge. Indeed, we shall see that the problem gets worse, not better, with 
this gauge choice. 

J jk 



Let us denote the polarization tensor for the GW in this gauge by tzfu- Then, so as to 



keep a basis with which we can compare with w^, we shall not change the definition of 

-~- -(si (s) 

the polarization vectors. Rather, we shall once again express Wj k in terms of by taking 

G$ ^ + ^r jk + 7 W W W , (Al) 

where 

t- - , ,(2L(2) / A n"i 

and aW, are coefficients which are to be determined. Clearly Wj k still satisfies the 

spatial components of the transverse condition. We still require that the two polarization 
tensors for the two polarization states of the GW be orthonormal to one another and this 
gives the constraint, 

1 = ( a (»))2 + + ( 7 W)2 , (A3) 



The 0-component of the transversality condition in eq. fli~7D can now be used with eq. (|A3|) 
to determine Wj k uniquely, 

Wj-fc = — sin ip cos i? cpjit — cos cos i9 + sin -d Wj k "* , (A4) 

where 

—t\ = q cos (f sin , 
-^-cc^-^ = g sin <p sin i? , 
Yp-w^ = qcosV, (A5) 
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and 

^(?H 2+ (?^) 2+ (?H 2 (Ae) 



For a RW universe, $ = 7r/2 while (p — 0. 

At this point, we make a few observations. First, the angles $ and ip are determined 
precisely by linear combinations of the terms found in eqs. ([54]) and (55). Second, these 



terms, which measures the infinitesimal change in the orthonormality conditions for the po- 
larization vectors, did not disappear when we change to this gauge choice. This underscores 
the fact that they are physically relevant. Third, the trace of the polarization tensors no 
longer vanish, 



E^ (+) , : = 2sin ^> E^ (X) i = 2 cos <^ cos??, (A7) 

3 j 



which will introduce additional terms to the GW lagrangian. 
In this gauge, the lagrangian for the GW now becomes 

-2R^ M h»h% - V p hV p h - hRTh^ , (A8) 

where h — hfc is the trace of h pv which does not vanish in this gauge and we are denoting 
the different choice of gauge with a hat. Once again we can divide this action up into three 
terms, the only difference being that the form of the gauge dependent piece is now different, 

h = \j V=gd 4 x (-V p hV p h - hR pu h pv ) . (A9) 

Proceeding in exactly the same way as before, we find that 

*k = \ E / ( E V P hsV p h s - E M S s,h s h s > - 2 E D^h^VhA , (A10) 

q V s ss' ss' / 



where M ss i and D pss i are defined in exactly the same way as M ss i and D pss i but with wfy 
replacing w^} in eq. fl£2]). 

It is straightforward to show that 
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V/V^r; = -2M 2 ++ - 2Dl x , 
V p u7 (x) ^Vr; = -2G = -4V M e« W 2) " . 

V^W^Vr; = -2H = -2 (V M eW - V M e<, 2) W 2)i/ ) . (All) 

Then 

M++ = M 2 + - V p <pV<p - sin2y?iJ + sin 2 y?/) 2 x , 
M 2 xx = M 2 ++ - V p Wtf - cos 2 #V p <pV p <p 
+ cos 2 d sin 2y?ff — sin 2d cos 

+ cos 2 $ cos 2 v^-D+x + (2 sin y^Vp^ — cos (p sin 2$V p (p)D p +x , 
M 2 x = sin Wy? - cos 2y? cos - sin i? sin y^G 

- (sin sin Wy? + cos y? cos $V p tf) D Hx + ^ sin 2<p cos $L> 2 x . (A12) 

While D p++ = D pxx =0 still, 

D p+X = sin^cos(/9-Dp +x + cos^V p y? . (A13) 

In a similar fashion, we find that 

T R = -\j2J2j a A d A xf ss ,h s ,h s , (A14) 

q ss' 

with 

T ++ = cos2v?T ++ , 

f x x = (sin 2 - cos 2 d cos 2y?) T++ , 

f +x = -cos^sin2v?T ++ , (A15) 

where we have used 

]T R\ m W^™T] = , 2 R^km^TT* = , 2 ^'kmOi = " T ++ ■ ( A1 6) 
ijfcrn ijfcm r/'fcm 

Finally, we have the gauge term, 
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2 

qss 



- E / a^xlw^w^V p h s V p h s , + 2m^V p z5 (s ' ) h s ,V p h s + 

qss' \ 

(^ ireH+ ? ( -^@4'l <ai7) 



where = cB^^/ a/2 and is given in eq. 

Combining the three terms together, we obtain 

i= \^l a " d ' x { - V P h s >V p h s 

qss' \ 

-2m (s) Vp£ {s % s ,V p h s - 2D pss ,h sl S7 p h s - { M ss , + f ss , + V ro (s) V p ro (s ' } 



^'^(S^M (A18) 

We can now see explicitly that instead of eliminating the coupling term, this choice of 
gauge merely makes things worse. In fact, the kinetic term for the two modes are now 
quite different than what we would expect and additional coupling terms now appear. It 
is, however, important to note that these additional terms have the same origin as those 
obtained by using the previous gauge. Namely, they all come from linear combinations of 
eqs. (54)- (|56|) and arise from the fact that the directions of the polarization vectors are 
continually changing with time. 
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